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HOMOGENIZATION AND CONTINUUM PERCOLATION 



DIMITRIS KONTOGIANNIS 



^—^ ' Abstract. We propose continuum percolation theory to study homogeniza- 

Cn ^ tion problems of elliptic equations. In particular, let g : R — > [0, 1] be a 

f^ I connectivity function that connects two points 3:i,X2 G K" with probability 

flj ' gil^l — 2:2 1), where [ ■ | denotes the Euclidean distance. Such functions have 

r^\ ' been introduced in continuum percolation theory. We involve the connectiv- 

ity function to study homogenization problems in random media. Our aim 
^ti' . is to improve and extend similar results that have been obtained for periodic 

domains. 



1. Introduction 



Homogenization theory is related to the asymptotic behavior of partial differen- 
tial equations describing physical phenomena in heterogeneous materials. In par- 
ticular, we look for the effective (homogenized) equations that describe the charac- 
teristics of the inhomogeneous medium as the length of a small parameter e tends 
CS| ' to zero. This small parameter is the length of the heterogeneity. 

^ I Some methods to approach such problems are, among others, the two scale 

expansion, the two scale convergence, the T— convergence. A summary of homog- 
enization techniques can be found in |10| . Relevant problems have been studied in 



m- 0, E], H, i, [13. 

In this note, we obtain averaged equations over randomly perforated domains. 
^-^ ■ For this purpose we use the connectivity function as defined in percolation theory 

[12]. More specifically, we consider the boundary value problem 



, - Au" -Xu" = f,xeG^{uj)nD 



?—( ' where the domain G^{uj) is a randomly perforated domain. We show that as the 

parameter e tends to zero, u^ converges to the solution of the problem 

, . Au- {X + c)u^ f,x e D 

^ ' u = 0,xedD 

The parameter c that appears in (1.2) is the extra term for the capacity which has 
been previously obtained for domains of periodic structure [7] as well as domains 
perforated by balls of random radius centered on the discrete Z" lattice [6]. Our 
approach is to provide a general model of randomly perforated domains so that (1.2) 
can be obtained without periodicity assumptions. The main feature of our model 
is the connectivity function which has been introduced in continuum percolation 
theory 12 . This function connects two points of a point process (Poisson process 
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for instance) with probability which depends on their distance. Thus, the main idea 
behind our model is 'as the density of points increases, points are connected only if 
their distance is of order e\ With this restriction, the model of random structures 
and shapes is based on the distance of the points. For example, we can think of thin 
structures such as tubules across the line segments or balls centered on the points 
with sufficiently small radius. To establish the homogenized equation, we use the 
ergodic theorem of stochastic integral functionals proved in [9] together with the 
ergodic properties of our model. 

We mention that in [T, the 'strange term' c appears in the equation when the 
'holes' have a critical size and examples for periodic domains are included. In some 
cases, the distribution c is given explicitly and the fundamental solution of the 
Laplacian on the annulus is an important tool. Our approach covers a wide range 
of non periodic perforated domains and critical values on the size of the structures 
are not under consideration in this work. It would be of interest though to find 
examples of random domains of this generality for which critical values can be 
approximated or computed. 

The paper is organized as follows: in section 2, we give a small review of results 
in stochastic homogenization, some of which we extend and improve in this note. In 
section 3 we include a brief review of discrete and continuum percolation theory and 
the construction of random models that we will use in homogenization problems. In 
section 4, we present our main homogenization results. Our work is an improvement 
of the homogenization method introduced by E.Khruslov(see for instance [H]). A 
key of the method is the definition of 'mesoscopic' characteristics, measureson a 
scale h with e « h « diam D and has been used in a variety of similar problems. 

2. Ergodic theory and stochastic homogenization 

Stochastic homogenization is based on the most generalized notion of periodicity, 
the stationarity. If, in addition, we have independence at large distances, we talk 
about ergodicity. Ergodic theory is related to the study of dynamical systems with 
an invariant measure. 

Let (f2, F, jjL, T) be a measure preserving dynamical system defined on a prob- 
ability space with the following structure: F is the a— algebra on fi, fi is the 
probability measure, T is the measure preserving transformation such that for any 
A e F, ^l{T-'^{A)) = ^(A). The system (rj,F,/Lt,T) is ergodic if the a-algebra 
of T— invariant events is trivial, that is, it occurs with probability zero or one. To 
see the importance of the theory on averaging problems, we state (among many 
versions) a subadditive ergodic theorem. 

A function // : A — )• R is called subadditive if for every finite and disjoint family 
(A,),67 with |A\U,e/A,| =0, 

^i{A)<Y,^i{A,) 

i 

We say that /i is dominated if < IJ.{A) < C\A\ for all sets A. Consider now the 
family of dominated, subadditive functions and the group of translations {tzIi){A = 
^(r^A), where t^A ^ {x &R" : x - z € A}. 

Theorem 1. (Ergodic) (see [1], [9]/' Let /i : 51 — > M" he a subadditive process, 
periodic in law, in the sense that ii{-) and TzfJ.{-) have the same distribution for 



HOMOGENIZATION AND CONTINUUM PERCOLATION 3 

every z e Z". Then, there exists measurable function (j> : il ^ R and a subset 
f2' C 51 of full measure such that 

lim — r^—. — = (h(uj) 
t-^oo \tQ\ ^^ ' 

exists a.e. lu E fl' and for every cube Q C M" . Furthermore, if /x is ergodic then (p 
is constant. 

We refer a paper from Dal Maso-Modica ^ for the use of ergodic theory in the 
calculus of variations, according to the following setting: 

Define the translation operator r^ that acts through the following relations: 
Tzu{x) = u(x — z), and t^A = {a; S M" : a; — z S A\, and the honiothety operator 
{peJ){u,A) = e'^J{peU,p^A) where {peu){x) = \u{£x), {peA) = {a; G K" : ea; G A}. 
A stochastic homogenization process is a family of random variables {J^)^^^ on a 
probability space {fl,F,P) that has the same distribution law with the random 
functionals given by [{p^J){uj)(u,A)] for u G W^'P{A). This means 

P{uj G n : F,{uj) G 5} == P{uj G Vl : PeJ{uj) G 5} 

for any open set S. 

We say that the random functional J is stochastically periodic, that is, J has 
the same distribution law as the random functional {tzJ){ijj){u, A) — J{tzU,TzA). 

Theorem 2. Let J(cj)('u, A) = / f{x,Vu)dx , where f satisfies standard growth 

J A 
conditions: k\p\'^ < f(x,p) < K\p\'^ for some positive constants k,K. Denote the 

minimizer of F by m{J, uq,A) — min{ J(w, A) : u ~ uq £ Wq ' (A)}. 

u 

If J is a random integral functional and if J and t^J — J{tzU,TzA) have the same 
distribution law, then the limit 

m{J{uj),uo,Qt) 



exists. If in addition J is ergodic, the limit is constant. 

The proof of this theorem is based on the ergodic theorem [T] , since the function 
p{uj){A) = m{J{u!),lp,A) is dominated and subadditive. 

Related works to stochastic homogenization include [1],[I1], [S]. We also men- 
tion the work of Cafarelli-Mellet [5^ in which the authors extended the results of 
Cioranescu-Murat [7] in the case that the obstacle problem is considered in a do- 
main perforated by balls of random radius centered at the Z" lattice. As they 
showed, depending on the capacity of the holes, we still have an additional term 
that appears in the averaged equations. 

3. Percolation theory and Random modelling 

Percolation theory deals with the behavior of connected elements in random 
graphs. A cluster is a simply connected group of elements. Percolation phenomena 
arise in transport, porous media, spread of deseases, conductivity problems, sea ice 
etc. The theory was introduced by Broadpent and Hammersley in 1957, when they 
considered the problem of fluid flow through a porous medium formed by channels, 
keeping in mind that some of the channels may be blocked. 
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In the discrete version of percolation, we consider the Z" lattice and for p e [0,1] 
we connect the point a: G Z" to each of its 2n nearest neighbors with probability 
p, independently of the other points. We describe distances in the lattice in the 
following sense: two vertices x, ?; G Z" are neighbors if |a; — y| = 1. We can also 
define boundary dQn of the sets Qn = [1, n]'' n Z", where dQn = {y E Z"^ : 3x ^ 
Qn-\x-y\ = 1}. 

According to this setting, each pair of neighbours has an edge between them 
with probability p. The edge is also called a bond. A path is a finite or infinite 
alternating sequence (21,61,2:2,62,..) of vertices Zi and bonds e^ such that Zi ^ Zj 
and Ci ^ Cj for i ^ j. Two vertices are connected if there is a finite open path from 
one to the other. An open cluster is a set of connected vertices that is maximal 
with respect to this property. It can be either finite or infinite. 

The main question that arises is if there exists a critical threshold of p at which 
an infinite cluster occurs. In many cases, like the two-dimensional lattices, the 
critical value Pc can be computed explicitly. This critical probability is an increasing 
function of p. 

4. Continuum Percolation 

A more general situation appears in the models of continuum percolation, where 
the integer lattice is replaced by a random set of points in M" . The random posi- 
tions are usually formed by the realization of a point process. Point processes are 
important models and have been used in a variety of problems such as environ- 
mental modeling, air pollution, weather radar images, traffic networks, statistical 
mechanics etc. For a complete account of the theory, please see [l2] . 

A point process is thought as a random set of points in the space. In particular 
we give the following definition: 

Definition 1. Let _B be a Borel subset of R" and denote by N the set of all 
counting measures in M". Let ■;/; e A^ be a counting measure, i.e. a measure which 
is 1 on each point x E B. Then N is identified as the set of all such configurations 
of points in M" without limit points. According to this setting, 'ip{A) — random 
number of points in A, for any set A C M". Then, a point process is defined as a 
measurable map X from a probability space (J7i, F, P) into {N, M, P), where M is 
the correspoding cr— algebra. 

The periodicity of the lattice-type structures is replaced with the assumption 
that the point process X{-) is stationary: 

If Ta is the translation in R" by a vector a, Ta{/3) = /3 -h a, V/3 G R", then T^ 
induces a transformation Sa '■ N —i' N through the operation {Saip){A) = iJ;(T~^A) 
VA E i?" and similar operation holds for the set-measures. 

Definition 2. The point process X is stationary if its distribution is S'^— invariant 
for any a E K". 

Two common models in the theory of continuum percolation are the Boolean 
model and the random connection model. Both models are based on occurences 
of Poisson processes. From this point, we assume that the point process X is a 
Poisson process of density A: 

(i) for any collection of mutually disjoint sets Ai,A2,...,Ak, 
the random variables X{Ai),X{A2), ...,X{Ak) are mutually independent, and 
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(ii) for any bounded set A C S" and any non- negative integer fc > 0, P{X{A) = 

k) = e^'*'''"*' -; , where l{-) is the Lebesgue measure. 

fc! 

4.1. The Boolean model. The Boolean model is driven by a Poisson process X 
and each point of the process is the center of a ball of random radius. The region 
in the space that is covered by at least one ball is called the occupied region and 
its complement is the vacant region. To be able to have shifting properties, we con- 
struct the models as follows. The Poisson process is defined in a probability space 
(r2i,Fi,Pi) and we consider a second space ^^2 = 1 [ 1 [ [0, oo) equipped with 

the probability measure /i on [0, cxd). Setting f2 = i7i x U,2 with product measure 
P — Pi X P2, the Boolean model is defined as the map {uji,lu2) — > {X{uji),ll!2) from 
fl into iV X il2. 

According to this construction, the radii of the balls are independent of the point 
process and we obtain the shifting properties that ergodic theory requires. 

We use the notation of [13] and we denote by (X, p, A) the Boolean model ob- 
tained from a Poisson process X of density A and radius random variable p. 

4.2. The Random connection model. As in Boolean models, the Poisson pro- 
cess is the first characteristic of the model and it assigns randomly points in the 
space. The second characteristic of the model is the connection function, which 
plays an essential role to the model and the homogenization process as we will see 
later. A connection function g : R+ -^ [0,1] connects two points xi,X2 £ X with 
probability g{\xi — X2\), where ] • | denotes the Euclidean distance. Depending on 
the construction that we need, we can choose g with specific characteristics. For 
example, we may assume that g is decreasing with respect to the distance. Such 
models are defined in product spaces as before. We denote them by {X,g, A). 

Definition 3. Two points x,y oi the process are connected if there is a sequence 
of points {a;o = x, a;i, X2, ...., Xn — y} such that each pair of points Xi, Xi+i are the 
endpoints of a line segment (edge) {xi, 0:^+1} for all i = 0,....,n — 1. As in the 
discrete percolation, a component is a set of points such that any two points of this 
set are connected and the set is maximal with respect to this property. 

4.3. Ergodic properties of point processes. The one dimensional ergodic the- 
orem is stated as follows (see for instance |12j): 

Theorem 3. Let (fl, F, fi,T) be a measure preserving dynamical system and let f 

_. n—l 

be p~integrable function on Vt. Then, — y /(T*(w)) -^ E{f\I){uj) , as n -^ 00 a.s. 

where I is the a— algebra ofT— invariant sets. 

For point processes, we identify any element lo £ fl with a counting measure. 
Then the shift by distance t, Tt, induces a transformation through (StUj){A) — 
uj{Tf^A) for all measurable sets A G K. 

Definition 4. A stationary point process is ergodic if the group {Sx ■ x £ M"} 
acts ergodically on {il, F, p). 

An important aspect of continuum percolation models is that the ergodic prop- 
erties of the process X are carried over the two models. The following results are 
well known [T2l: 
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Theorem 4. ; A Poisson point process is stationary ergodic. 

Theorem 5. Suppose that the point process X is ergodic. Then the random con- 
nection model {X,g,X) and the Boolean model (X, p, A) are also ergodic. 

Let (r^i , Fi , Pi) be a probability space associated with the Poisson process X and 
let cj G ill be a point configuration on R" that is assumed to be locally finite and 
countable. This means that we have a finite number of points hits every compact 
set K C M" almost surely: 

P{uj e n : tj){K) < cx) for all compact K C M") = 1 

5. Construction of random domains for homogenization 

5.1. Random connection models. We start using the random connection model 
{X, g, A) in the following way: 

Suppose that a; is a given realization for X which is locally finite and let Xi E X 
be a given point of this realization. 

Consider the annulus A = [x (z K." : ci < |a; — x^j < C2}, where ci,C2 are positive 
constants with ci < C2. 

We want to connect the point Xi with all the points in A that are given from X. 
For this purpose we choose the connection function 

J 1 if Ci < \X - Xi\ < C2 

g(\x — Xi\) = < 

I otherwise 

For a point Xj G A, we denote by hji^) ~ l{xi,Xj) the line segment with 
endpoints Xi,Xj and let Tc-^/2{hj){'^) the tube of radius ci/2 surrounding lij. Let 
now T{xi){uj) = UjTcj/2(^ij)('^) ^^^ F(uj,ci/2) = UiT(xi){uj) for all points Xi of 
the process. 

Thus, the set F{uj, ci/2) is the union of random tubes obtained from the given 
realization of the point process. Let G(a;,ci/2) = M" \ F{uj,ci/2). We define the 
indicator function 

a(oj, a;) = 1 - min{Xp(^^^ci/2)A} 

which is zero in the union of tubes and one elsewhere. 

Let D be an open, bounded domain of R" and consider the random functional 

J{uj){u,D) = / a{uj,x)\Vu\^dx ^ f \\7u\^dx for u e W^-^{D). This func- 

Jd JG(uj)nD 

tional is periodic in law and independent at large distances, thus ergodic. 

Furthermore let {p^J){u^A) — e''^J{p^u,p^A) where {p^u){x) — -u{ex), (peA) — 

e 
{x e R" : ea; G A}. Then the family 

J%u,D)^ p,J{u,D)^ I iVwpdx 

JE"G{ui)nD 

satisfies the assumptions of theorem 2. Note that the pe— homothetic functional is 
the functional obtained if we scale by e the distance between the connected points 
of the set F{uj,ci/2) that corresponds to the union of tubes eF = F{e(jj,eci/2), 
where euj maps to the point measure whose support is {exi} and {xi} is the support 
of X{uj). Note that the scaling properties of this model are the same (in terms of 
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distribution) with the model that we have if we choose 

J 1 if ci£ < la; - Xi\ < C2e 
ge(\x — Xi\) = < 

I otherwise 

with density function X/e. Let us define F^{uj) = sF = F{euj,eci/2) and G^(w) = 
eG{ijj) — K." \ F'^{ijj). According to this model, only points of £— distance are 
connected. 

Remark 1. The radius of the tubes need not necessarily be constant. We may 
consider the tubes Tp^(^^^'^{l{xi,Xj)) = {x G 17 : d{x,l(xi,Xj)) < pe{x,uj)}, where 
Ps{x,u;) is continuous function on x, stationary ergodic and for some positive con- 
stants pi.e,P2.e, Pi.e < PeixjUj) < p2.e- This is duc to the fact that the product of 
two ergodic processes is also ergodic. 

5.2. Boolean models. We take again a probability space {ili,Fi, Pi) associated 
with the Poisson process X and let a; £ O be a given point configuration on R". 
We consider a second space ^2 = TT TT [0, 00) equipped with the probability 

measure p on [0, cxd) for the sequence of independent, identically distributed random 
radii. Let us denote by w = (w, f) the realization of this Boolean model in O = 
r^i X fl2, where f = (ri, r2, ...). 

Let F{lo) be the union of random spheres obtained from the given realization of 
the point process. Let G(oj) = K" \ F{uj). We define the indicator function 

a(cij,x) = 1 — min{Xp (^1^-^,1} which is zero in the union of spheres and one 
elsewhere. 

Let D be an open, bounded domain of R" and consider the random functional 

J{uj){u,D)— / a(a;,x)|VM| dx — jVul dx. This functional is periodic in 

Jd JG{uj)nD 

law and independent at large distances, thus ergodic. 

Then, the family p^J{uj){u^D) ~ J^{u,D) satisfies the assumptions of theorem 
2. 

A Boolean model gives, in general, a union of spheres which may intersect. One 
way to model non-intersecting spheres is to combine the Boolean and the random 
connection model in the following way: 

Suppose the RCM (X, 5,A) is applied on a bounded region of M". We want to 
assign every endpoint of the line process as the center of a ball of random radius. 
For fixed e > {) there is a set of points from the point process X. In our case, 
instead of constructing tubes, we let every point be the center of a ball with radius 
p{uj) < rmnd{xi,Xj){uj), where the minimum is taken over all the pairs of points 
X oi X{llj). Note that, without any affect to our proofs, we may assume that p{uj) 
is identically distributed random variable taking maximum value rmnd{xi,Xj){u;). 
We consider for simplicity the first case. According to this construction, we obtain a 
domain randomly perforated with balls of radius and with positive minimal distance 
to each other. Let us define 



F%Lo)^ \jB{ep{Lu),ex,)nD 



i>i 



and G^ (to) n D = D \ F^ (co). Note that mesi^''(a;) tends to zero as e ^ 0. 
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6. The Dirichlet problem 



6.1. Boolean models and domains with fine-grained boundary. Let w G Oi 

be such that X{uj) is locally finite. 

We consider the Dirichlet problem of the form 

, . Au^ -Xu" = f,xeG%Lu)nD 

^ ' u' ^0,xedG'{Lu) 

for u' e W^'^{G^{u})), f e L'^{D). Here, G^{uj) is the domain perforated by 
non-intersecting balls. 

To define capacity characteristics for the massiveness of Bf{uj), let us consider 
the quantity 



cap(B) = inf / \Vv\'^dx 



for closed, bounded sets B in R", over all functions v G C^{W^) taking value 1 in 
B. We consider only the case n > 3 so that cap(i3) defines the Newton capacity. 
Note that the capacity is invariant with respect to translations and rotations. In 
addition, if B^ is the e— homothetic contraction of B, cap(i3e) = e""^cap(_B). 

We also note that, clearly, as e ^^ 0, the diameter of the balls tend to zero. Now, 
note that the limits lim \^ cap(i?f (w)) = cap(G'^(a;)) = C exists due to the ergodic 

[D) 

properties of the model and theorem 1. The sum is taken over all balls strongly 
contained in D. 

We state the homogenization theorem for the Dirichlet problem related to Boolean 
models. 

Theorem 6. The family of solutions u'^ of the Dirichlet problem (3) (extended 

n{e) 

by zero in M Bl{u))) converges in L^{D) to the solution u of the boundary value 

i=l 

problem 

, , Au~{X + C)u^ f,x€ D 

^^ > u^O,x€dD 

The proof of this theorem can be found in [llj under general assumptions which 
have been proved in our cases. 

6.2. RCM and connected domains. Let -F'^(w) be the random set of channels 
that we constructed in section 4.2. We assume that the Poisson process X is locally 
finite in the sense that a finite number of points hits every compact set K C M" 
almost surely. That is, P{uj G J7 : ip{K) < cxd for all compact K C M") = 1. 

Let G'^(w) — D\F^{uj). We want to show that for sufRciently small e, the volume 
of G'^{uj) is strictly positive with probability 1. That is, there is e{uj) > such that 

P{lj en-. \G%uj)\ > for all e < e{u)) = 1 

For simplicity and without loss of generality, suppose that n — 2 and D is a 
square of size 1/e > 0. Consider a partition of D into squares Di, i — 1, .., 1/e^ of 
size 1. Suppose that the point process X is applied in D. Then, the probabilty of 
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having zero points in a given Di is e^^ > 0. Aceording to the law of large numbers, 
if li{uj) is the indicator function of the empty square Di of the partition, 



1 

P{Loen: lim -— - y U(uj) = e-^) = 1 

e^oo l/£2 Z^ ^ ' 
' i—1 

which means that, in the limit, there are empty squares with probability one, as 
needed. Since the points, within an e— range, are connected, the set 

{^.3} {kd}^{t.j} 

(the intersection of the tubes) is non-empty, K^{uj) C F^{u!). Clearly, as e -^ 0, 
mesF^(w) — ?■ and G^{uj) becomes denser in il. One may consider two different 
homogenization problems: the first is when the elliptic equation is defined in ^."^(a;) 
for M^ G W^''^{G^{u!)) and the second is the problem of decreasing volume when 

We consider the Dirichlet problem of the form 

..ON -Au' + Xu' = f,xeG%uj) 

^ ' u" ^0,xedG%Lu) 

for w^ G VF^'^(G'^(a;)), / G L'^{D) and as £ -> 0, G"^(cj) is approximately D. 
In general, F^{lu) consists of connected components, but it is not necessarily a 
connected set. 

Standard elliptic theory gives the existence of solutions. We extend the functions 
u'^ by zero in F'^(uj) and we keep the same notation for the extended sequence of 
functions. We denote by Q% the n-cube centered at x of length h, diamD >> h » 
e, and we define the local capacity functional 



(6.4) cap(x, /i,£,a;) =inf / \\7u''{y)\^dy 

over aU u^ G W^'^{Ql) : m-^ = in F^lu). Clearly, if F^i(cj) nQ^C F^^lo) n Ql, 
then cap(x, /i, £i, cj) < cap(x, /i,£2, w). Thus, the capacity functional measures the 
massiveness of F'^{uj) in il. Note that, from our construction and due to theorem 
(5), F'^{u;) and its complement are periodic in law and ergodic in the sense that 
disjoint cubes have the same distribution. 

Note that cap(x, /i, £, cj) = inf / a{x,uj)\'Vv^\'^dx, where 

a{x,uj) = 1 -min{l,XTp(;.^)} 

which is in the union of random channels with endpoints the points of the Poisson 

process and 1 elsewhere. 

Hence, theorem 2 is applicable over the class of functions in W^'^{D). Thus, the 

constant limit 

ca.p(x,h,e,uj) 
c — lim lim 

h^Oe^O ft," 

exists almost all a; G ili. Hence, we can also assume that 

ca.p{x, h,e,uj) 
lim sup < A 

for all a; G Z? with A independent of h. Our main homogenization theorem is the 
following: 
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Theorem 7. Let F^{u)) he the sequence of RCM domains constructed in section 5 
and G^{u]) he its complement set. Let u^ G W^''^{G^{uj)) he the family of solutions 
of the houndary value prohlems (6.3) extended by zero in F^(uj). Then, as e — > 0, 
u^ converges in L^{D) to the limit u G W^'^{D) which solves the boundary value 
problem 

, , Am - (A + c)u = f,x e D 

^ ' u = 0,x£dD 

Proof. Note that the (extended by zero) solution u'^ of (6.3) is the miniinizer in D 
ofthc functional n[M^] = / \Vu^\^ + X\u^\^ +2fu^dx = f |Vm^P + A|m^P + 

JG'(Lo)r\D Jd 

2fu^dx over the class of functions u'^ e W^''^{G^{uj)). 
Thus, T^[u^] < T^[0] = which implies that 

\Vu'\^ + X\u'\^dx <2\\u'\\L2^n)\\f\\mn) 

G^{uj)nD 

Using the Friedrich's inequality ||w'^||i2(£)-) < C|| VM'^||L2(i5) we obtain that 

where C is independent of e. Since u^ in bounded, it has a subsequence, still de- 
noted by M*^, that converges weakly in W^'^{D) and strongly in L^{D) to some 
function u £ W^'^{D). 
To prove the theorem, it is enough to show that the limit u is the minimizer of the 

functional T[u] = / |V7ip + (A + c)\u\'^ + 2fudx. 
Jd 

Stepl: We first estabhsh the inequality hm sup F"^ [u*^] < r[w] for all w G W^'^{D). 

For this purpose, we consider a partition of D with cubes Q" ~ Q{x" , h) centered 
at x" of size h, so that UaQix", h) is a cover of D and the points x" form a periodic 
lattice of period ft, — r, r to be chosen. A partition of unity {(pa} of C^ functions 
subordinated to this covering is given as follows 

(1) <(/.„< 1 

(2) </.„ = if a; ^ g", ./)„ - 1 if X e Q" \ Up^^Q^ 

(3) ^0a(a;) = 1, if xeD 

(4) \V4,a\<C/r 

Let us denote by w" = ti"('^' the minimizer of cap(x, h, e, cj) in the cube centered at 
x". For w G C'^{D), compactly supported in D, define 

n{h) n(h) 

(6.6) wl{x) = Y, w(x)w"(x)0„(x) ^ w{x) + Y, w{x)[v°'{x) - l]</-„(x) 

Q— 1 a— 1 

so that w% e W^^'^iD) and w%{x) == in F^iu). Thus, F^[m'^] < T^[w%]. Under our 
assumptions, 

(6.7) / |Vw"pdx<C/i" 
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We denote by Qh ~ Qh \ ^fi^aQ^ the concentric cube centered at x" of size 
h^h-2r. Then, 

( \Vv°'\^dx 

= f iVv^'l^dx- f \Vv''\^dx + 0{rh"-'^) 
< ca.p{x,h,e,uj) — cap(a;, /i, e, w) + 0{rh"'^^) 

Choosing r = h^^''''^ = o(h), where < 7 < 2 is a positive parameter, we obtain 



(6.8) / |Vt;"|^dx = o(/i") 

Note that Fricdrich's inequaUty gives 



(6.9) / |w" - l^da; < C(r) / {Vv^l^dx = o{h''+') 

Differentiating (5.6), we have 

a a a 

We substitute (6.10) into r^[ii;|] to obtain 

T-[wl\^ [ \Vwi\^ + X\wl\^ + 2fwldx^ f \Vw\' + Xv\w\' + 

n(h) 5 

+ 2fwdx +J2 |Vw"pw2,/,2 da; + J2 U{,e, h). 

a=l''Qh i=l 

where 



N(h) 

2 

'Qh 

L2{e,h) = 



Li(e,/i)== V 2 / {f + \w){v°' ^ l)w(j)o,dx 

Jq?: 



Z^ JQh^qI it ^^^' " dx.'^jydx, ^ dx,'^ 



+ Aw2 0„0^}(w" - l)(w'3 - l)da; 



Ar(/i) „ 



C.R ^=l JQtr^Qi ^^» ^^» 



L.{e,h) =f^±2 j ^u,<^^'M-^dx 
^-^ ^-^ /n° oxi axi 
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Taking into account the properties of <j)a, (6.7) — (6.9) and the fact that the 
number of cubes Qf^ which intersect Q^ is no more than 3" , we have 



lim Um > Li(e, h) — 

i=l 

From the properties of smooth functions (pa we now have 

w'^(j)a\Vv°-\^dx < CiD^cap(a;",/i,e,w) 



where Wa is the mean of Wa over Qf^. Summing over the cubes and letting s tend 
to zero we have 

(6.11) hmSUp^^o Ea /q^ I V-y" 12^202 gf^. < j^ cw2d2; + 0{h) 

Combine these inequahties to see that hmsupr^[u^] < T[w] for all twice differ- 

entiable functions with compact support in D. Using a density argument, this 
inequality holds for all w € Hq{D). 

Step 2: To show the reverse inequality liminf r'^[u^] > T[u\, pick a sequence 

us{x) G Cq{D) such that \\us — u\\[{i/ij\ < e, where u is the weak limit of the 
sequence of minimizers u^ of r'^[-] in Hq{D). 

According to lemma 3.2 on JTj pg.73, there is a sequence {uf} G Hq{D,F'^{uj)) — 
{v e Hq{D) : t) = in F^(w)} that converges to us and satisfies \\ul — u"^ \\ h^ (d) < 

C\\us -u||ffi(D). 

Take now the cubes Q^ that belong to the set Dg = {x ^ D : \us{x)\ > S} for 

y£ 

positive parameter S and in each of these cubes define the function i;^ = — so that 
vl -^ 1 weakly in L2(Qa). 
Clearly, 

(6.12) /q„ \Vv%\'^ + h-^-^\vl^-l\'^dx>c&^{x'',e,h,-i,uj) 
and 

lo io\ 2< _ J_*4 )_du5_ _ "j-fj dus 

V / dxi Us dxi Us dxi ul dxi 

Using the expansion (6.13) in (6.12) and taking into account that Ug — > us strongly, 

we get 

(6.14) 

r ■ I |l2 

/q° |Vu||^dx>cap(a;",e,/i,w)[minQ^ \us\? + ^^-^^ Jqc, \Vus\^dx - G{e,S,S,h) 
where lim G{e, S, S,h) — for fixed h, S, S. 

e— >0 

Finally, we sum over all cubes that intersect G'^(w), 

nK]>f: / |Vu.|2 + f] ^^P("".f'^'") [supK|]2fe"+ / AK|2 + 2M 

" [supq„ |u5|]2 - [miuQo |m5|]2 f- 
-> ^ ^ / \Vu,\^dx-NG{e,5,5,h) 
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We let h ^ for fixed S to obtain 

limr=[u^]> / \Vus\'^ + culdx + / \uj + 2fu'jdx 

^-*° JDg Jd 

Let now 6, S, e tend to zero: T[u] < T[v] for all v £ Hq{D). D 

7. The Neumann problem, strong and weak connectivity 

In [llj . the Neumann problem for elliptic equations is considered. The conditions 
of convergence are formulated in terms of mean local characteristics, that we will 
describe. 

Let D C M" be a bounded domain and G^{uj) fl D be the domain constructed 
in section 5. We seek the compactness of the sequence u'^{x) defined in G^{Ld) n D, 
i.e, the possibility that a subsequence u^''{x) converging to some u £ L^iD) in the 
sense that 

Wu"" -u||L2(G.(^)n£i) -^ 
as e ^ 0. 

The domains {G"^(aj)}e>o are said to be strongly connected if any sequence of 
function u^ defined in G^(aj) is compact with respect to the last convergence. 

The main local characteristic of the medium is the functional 

PlhiO^ ,„, inf ^ ^ ^^ f{\Vv-\^ + h-'-^\v-^ix~z,0\^}dx 
Note that the function u| which minimizes P^ii{£,) for any ^ G M" can be written 

n 

in the form vf = y £,ivf, where vf is the corresponding minimizer for ^^ — ei. 

1=1 
Thus, if we write 

aij {z, £, h,uj) ^ / Vwf • Vu| 

+ h-'-'ivt - (x. - COM - i^j - ^j)]dx 

n 

then P^ hiC) takes the quadratic form P^ ii{^) = 2. o,ij{z,e, h,u})$,iS_j. 

The matrix [a.y] is the local mean conductivity tensor of the medium at the point 
z. 

Similar capacity-type functionals can be defined for weakly connected domains. 
Such domains consist of a finite number of strongly connected components. We 
note that the domains G^{uj), e > are weakly connected domains, since the 
connectivity function g{\ ■ |) gives, in general, only connected components. 

Furthermore, the temperature distribution w^ of the porous medium is the min- 
imizer of 

JG'{u) 

over the class {u" 6 VK^'^(G^(w)) : u^ = uq on d{D)}. 

Suppose for the moment that u'^ can be extended to vf such that 

||'u^||h'i,2(d) < C uniformly with respect to e (see [15j for improvements). Then, 

up to a subsequence v!^ converges to a function u in L'^{D), which is almost linear 
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to every sufficiently small cube K^, i.e. u{x) = u{z) + {x — z, Vu{z)) + o{h?). Then 
for e small enough, 



/ 

JQ: 



\v'{x) - (x - z, Vu{z))\^dx = 0(/i"+^) 
/QjnG-(c^) 

where v^{x) — u^{x) — u{z). Thus, we can assume that the function v'^(x) is the 
minimizer of 

Jl{v')=j \Vv^^dx = 0{h-) 

Under this consideration, 

J'nivn ^ Plhi^uiz)) ^ o{h-) 

as /i — ^ 0. 

The construction of G"^(w) says that J^ is periodic in law and independent at 
large distances in the following sense: 

Define the function a{uj,z) — 1 — min{l,X7i (h)} which is in the union of 
random channels with endpoints the points of the Poisson process and 1 elsewhere. 

Consider the random functional J{uj){u, A) ~ / a{uj,x)\'Vu\ dx — I |Vm[ dx 

J A Jg{ui) 

and define p^J — J^(w). 

Thus, theorem 2 shows that lim , ^ , exists. This implies that 

h^o \Ql\ 

a^j{z,£,h,uj) _ 

hm hm — —z = aij(x) = a 

exists almost all w € fii and it is constant. 

It is worth to mention that, in order to apply the same homogenization technique, 
a major open problem is to prove the extension of functions from randomly perfo- 
rated domains for Neumann problems, since the strong connectivity assumption of 
[11] does not hold with probability one. 

8. Homogenization for sets of decreasing volume 

In [11], [13] the homogenization problem for sets of decreasing volume is con- 
sidered. In such problems, we allow the measure of the domains of definition of 
functions to vanish. A main assumption to define appropriate type of convergence 
is a uniform density of the vanishing set within the domain D: li F^ <Z D \s such 
that |_F^| — J> as e — !> 0, then for sufficiently small e > we require 

Cr"|i^"| < \B{x,r)r\F''\ < C-V"|i^^| 

for every ball B{x, r) centered at x G 13 of radius r > 0. 

In particular, let us consider the set F^{oj) of tubes as in section 5, so that 
1^*^(0;) I — 7- as £ — > 0. We want to show that the density assumption holds for 
these domains with probability 1. Recall that F'^{uj) — eFiuj) is the union of tubes 
with endpoints the points of X which has the same distribution with the model of 
density A/e and height at least cie and at most C2S. Thus, as e decreases, the mean 
density increases and the distance between the pairs of connected points decreases. 

Suppose that the last density condition does not hold. This means that there is 
a ball B{x, r) such that for all e < e(w) and w G ri2 C fii with \VL2\ ^ 0, 

\B{x,r)r\F^{oj)\ = 
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This means that, for aU small e, there are no points Xi,Xj G B{x,r) of the process 
X with ci£ < \xi ~ Xj\ < C2£. On the other hand, 

P{X{B{x,r)) = 0) = e-^''"/'^ ^ 

as e — !• whieh is a contradietion. 

9. Further remarks 

The modeling based on connectivity functions may be used for more general 
shapes and random structures, since the total measure of the sets that it produces 
is forced to tend to zero. 
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